Search for: All records

Isospectrality is a general fundamental concept often involving whether various operators can have identical spectra, i.e., the same set of eigenvalues. In the context of the Laplacian operator, the famous question “Can one hear the shape of a drum?” concerns whether different shaped drums can have the same vibrational modes. The isospectrality of a lattice in $d$-dimensional Euclidean space ${\mathbb{R}}^d$is a tantamount to whether it is uniquely determined by its theta series, i.e., the radial distribution function $g_2\left(r\right)$. While much is known about the isospectrality of Bravais lattices across dimensions, little is known about this question of more general crystal (periodic) structures with an $n$-particle basis ( $n\ge 2$). Here, we ask what is $n_{\text{min}}\left(d\right)$, the minimum value of $n$for inequivalent (i.e., unrelated by isometric symmetries) crystals with the same theta function in space dimension $d$? To answer these questions, we use rigorous methods as well as a precise numerical algorithm that enables us to determine the minimum multiparticle basis of inequivalent isospectral crystals. Our algorithm identifies isospectral four-, three- and two-particle bases in one, two, and three spatial dimensions, respectively. For many of these isospectral crystals, we rigorously show that they indeed possess identically the same $g_2\left(r\right)$'s for all values of $r$. Based on our analyses, we conjecture that $n_{\text{min}}\left(d\right)=4$, 3, 2 for $d=1$, 2, 3, respectively. The identification of isospectral crystals enables one to study the degeneracy of the ground-state under the action of isotropic pair potentials. Indeed, using inverse statistical-mechanical techniques, we find an isotropic pair potential whose low-temperature configurations in two dimensions obtained via simulated annealing can lead to both of two isospectral crystal structures with $n=3$, the proportion of which can be controlled by the cooling rate. Our findings provide general insights into the structural and ground-state degeneracies of crystal structures as determined by radial pair information. Published by the American Physical Society2024 

The probability of finding a spherical “hole” of a given radius r contains crucial structural information about many-body systems. Such hole statistics, including the void conditional nearest-neighbor probability functions GV(r), have been well studied for hard-sphere fluids in d-dimensional Euclidean space Rd. However, little is known about these functions for hard-sphere crystals for values of r beyond the hard-sphere diameter, as large holes are extremely rare in crystal phases. To overcome these computational challenges, we introduce a biased-sampling scheme that accurately determines hole statistics for equilibrium hard spheres on ranges of r that far extend those that could be previously explored. We discover that GV(r) in crystal and hexatic states exhibits oscillations whose amplitudes increase rapidly with the packing fraction, which stands in contrast to GV(r) that monotonically increases with r for fluid states. The oscillations in GV(r) for 2D crystals are strongly correlated with the local orientational order metric in the vicinity of the holes, and variations in GV(r) for 3D states indicate a transition between tetrahedral and octahedral holes, demonstrating the power of GV(r) as a probe of local coordination geometry. To further study the statistics of interparticle spacing in hard-sphere systems, we compute the local packing fraction distribution f(ϕl) of Delaunay cells and find that, for d ≤ 3, the excess kurtosis of f(ϕl) switches sign at a certain transitional global packing fraction. Our accurate methods to access hole statistics in hard-sphere crystals at the challenging intermediate length scales reported here can be applied to understand the important problem of solvation and hydrophobicity in water at such length scales. 

The knowledge of exact analytical functional forms for the pair correlation function g2(r) and its corresponding structure factor S(k) of disordered many-particle systems is limited. For fundamental and practical reasons, it is highly desirable to add to the existing database of analytical functional forms for such pair statistics. Here, we design a plethora of such pair functions in direct and Fourier spaces across the first three Euclidean space dimensions that are realizable by diverse many-particle systems with varying degrees of correlated disorder across length scales, spanning a wide spectrum of hyperuniform, typical nonhyperuniform, and antihyperuniform ones. This is accomplished by utilizing an efficient inverse algorithm that determines equilibrium states with up to pair interactions at positive temperatures that precisely match targeted forms for both g2(r) and S(k). Among other results, we realize an example with the strongest hyperuniform property among known positive-temperature equilibrium states, critical-point systems (implying unusual 1D systems with phase transitions) that are not in the Ising universality class, systems that attain self-similar pair statistics under Fourier transformation, and an experimentally feasible polymer model. We show that our pair functions enable one to achieve many-particle systems with a wide range of translational order and self-diffusion coefficients D, which are inversely related to one another. One can design other realizable pair statistics via linear combinations of our functions or by applying our inverse procedure to other desirable functional forms. Our approach facilitates the inverse design of materials with desirable physical and chemical properties by tuning their pair statistics.